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RECENT PROGRESS IN LEPTOGENESIS* 

W. BUCHMULLER 
Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany 

f— ^ ' After recalling the general virtues of leptogenesis we compare two realizations, 

£""v ' Affleck-Dine leptogenesis and thermal leptogenesis, which generically lead to dif- 

^«. 1 t ferent predictions for neutrino masses. Finally, we describe the progress towards a 

full quantum mechanical description of the basic non-equilibrium process of lepto- 
genesis beyond the approximations involved in Boltzmann's equations. 
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C*~) \ 1 Why Leptogenesis ? 



One of the main successes of the standard early-universe cosmology is the prediction of the 
abundances of the light elements, D, 3 He, 4 He and 7 Li. Agreement between theory and 

#-/-\ ' observation is obtained for a certain range of the parameter r\, the ratio of baryon density 

l/~) • and photon density^, 
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where the present number density of photons is n 7 ~ 400/cm 3 . Since no significant amount 
of antimatter is observed in the universe, the baryon density yields directly the cosmological 
baryon asymmetry, Yg = (ns — n§)/ s ~ 77/7, where s is the entropy density. 

A matter- antimatter asymmetry can be dynamically generated in an expanding uniyerse 

Qh, if the particle interactions and the cosmological evolution satisfy Sakharov's conditionsEl,i.e., 

CD 
i-C • baryon number violation 



C and CP violation 



• deviation from thermal equilibrium 
C3 ' 



Although the baryon asymmetry is just a single number, it provides an important relation- 
ship between the standard model of cosmology, i.e., the expanding universe with Robertson- 
Walker metric, and the standard model of particle physics as well as its extensions. 

At present there exist a number of viable scenarios for baryogenesis. They can be classi- 
fied according to the different ways in which Sakharov's conditions are realized. Already in 
the standard model C and CP are not conserved. Also baryon number (B) and lepton num- 
ber (L) are violated by instanton processes!. In grand unified theories B and L are broken 
by the interactions of gauge bosons and leptoquarks. This is the basis of the classical GUT 
baryogenesis! Analogously, the L-violating decays of heavy Majorana neutrinos lead to 
leptogenesisB. The initial abundance of the heavy neutrinos may be generated either ther- 
mally or from inflaton decaysB'Q. In supersymmetric theories the existence of approximately 
flat directions in the scalar potential allows for new possibilities. Coherent oscillations of 
scalar fields may then generate large asymmetries^. 
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The crucial departure from thermal equilibrium can also be realized in several ways. 
One possibility is a sufficiently strong first-order electroweak phase transition^. In this 
case CP violating interactions of the standard model or its supersymmetric extension could 
in principle generate the observed baryon asymmetry. However, due to the rather large 
lower bound on the Higgs boson mass of about 115 GeV, which is imposed by the LEP 
experiments, this interesting possibility is now restricted to a very small range of parameters 
in the supersymmetric standard model. In the case of the Affleck-Dine scenario the baryon 
asymmetry is generated at the end of an inflationary period as a coherent effect of-scalar 
fields which leads to an asymmetry between quarks and antiquarks after reheatinglilJ. For 
the classical GUT baryogenesis and for leptogenesis the departure from thermal equilibrium 
is due to the deviation of the number density of the decaying heavy particles from the 
equilibrium number density. How strong this deviation from thermal equilibrium is depends 
on the lifetime of the decaying heavy particles, and the cosmological evolution. Further 
scenarios for baryogenesis are described in refJli. 

The theory of baryogenesis involves non-perturbative aspects of quantum field theory 
and also non-equilibrium statistical field theory, in particular the theory of phase transitions 
and kinetic theory. A crucial ingredient is the connection between baryon number and lepton 
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Figure 1. One of the 12-fermion processes which are in thermal equilibrium in the high-temperature phase 
of the standard model. 



number in the high-temperature, symmetric phase of the standard model. Due to the chiral 
nature of the weak interactions B and L are not conserved. At zero temperature this has no 
observable effect due to the smallness of the weak coupling. However, as the temperature 
approaches the critical temperature TE\^-f>i the electroweak transition, B and L violating 
processes come into thermal equilibriun£3. 

The rate of these processes is related to the free energy of sphaleron-type field config- 
urations which carry topological charge. In the standard model they lead to an effective 
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Figure 2. Sphaleron rate in Bodeker's effective theory, two lattice-Hmplementations of the HTL effective 
theory, and pure lattice theory interpreted as HTL effective theoryej. 

interaction of all left-handed quarks and leptonstj (cf. fig. [y), 

O b +l = JJ {qLiquqLihi) , (2) 

i 

which violates baryon and lepton number by three units, 

AB = AL = 3 . (3) 

The evaluation of the sphaleron rate in the symmetric high-temperature phase is a 
complicated problem. A clear physical picture has been obtained in Bodeker's effective 
theorytJ according to which low-frequency gauge field fluctuations satisfy the equation of 
motion 

D x B = o-E - C ■ (4) 

Here £ is Gaussian noise, i.e., a random vector field with variance 

(C^xX^x')) = 2*6 ij 6(t-f)5(yL-x') , (5) 

and a is a non-abelian conductivity. The sphaleron rate can then be written aaiil, 

r ~ (14.0 ± 0.3) -(o^T) 5 . (6) 

a 

A comparison with two lattice simulations is shown in fig. 0. From this one derives that B- 

and L-violating processes are in thermal equilibrium for temperatures in the range 

T EW - 100 GcV < T < T SPH - 10 12 GcV . (7) 



Sphaleron processes have a profound effect on the generation of the cosmological baryon 
asymmetry, in particular in connection with the dominant lepton number violating interac- 
tions between lepton and Higgs fields, 

Cal=2 = \h iL<t>ci Lj <t> + h.c. . (8) 

Such an interaction arises in particular from the exchange of heavy Majorana neutrinos. In 
the Higgs phase of the standard model, where the Higgs field acquires a vacuum expectation 
value, it gives rise to Majorana masses of the light neutrinos v e , v^ and v r . 

One may be tempted to conclude from eq. (Q) that any B + L asymmetry generated 
before the electroweak phase transition, i.e., at temperatures T > Tew, will be washed 
out. However, since only left-handed fields couple to sphalerons, a non-zero value of B + 
L can persist in the high-temperature, symmetric phase in case of a non-vanishing B — 
L asymmetry. An analysis of the chemical potentials of all particle species in the high- 
temperature phase yields a relation between the baryon asymmetry Yg ^Jn B — n B )/ s and 
the corresponding B — L and L asymmetries Yb-l and Yj,, respectively^, 

Y B = a Y B _ L = -i- Y L . (9) 

a — 1 

The number a depends on the other processes which are in thermal equilibrium. If these are 

all standard model interactions one has a = 28/79. If instead of the Yukawa interactions of 

the right-handed electron the AL = 2 interactions (ra) are in equilibrium one finds a = —2/3 

From eq. (f|) one concludes that the cosmological baryon asymmetry, if generated before 
the electroweak transition, requires also a lepton asymmetry, and therefore lepton number 
violation. This leads to an intriguing interplay between Majorana neutrinos masses, which 
are generated by the lepton- Higgs interactions (|8|), and the baryon asymmetry: lepton 
number violating interactions are needed in order to generate a baryon asymmetry; however, 
they have to be sufficiently weak, so that they fall out of thermal equilibrium at the right 
time and a generated asymmetry can survive until today. 

The connection between baryon and lepton number at high temperatures leads to very 
attractive new mechanisms to generate the cosmological baryon asymmetry. This includes 
large lepton number violating classical fields generated during an inflationary phase and the 
out-of-cquilibrium decays of heavy Majorana neutrinos. In the following we shall compare 
these two versions of leptogenesis which generically make different predictions for neutrino 
masses. 

2 Affleck-Dine leptogenesis 

In supersymmetric theories the AL = 2 lepton- Higgs interactions (|8|) are contained in the 
superpotential 

W - ^—^L^U , (10) 

ZMi 

where Hi and Li, i — 1...3 denote Higgs and lepton superfields in a particular basis, 
respectively. The supersymmetric standard model possesses the D-flat directionEZ 
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Figure 3. Contour plot for the baryon asymmetry Yb 
Tfl[GeV](vertical)-plane. The full lines correspond to Yb - 
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ub/s in the rriv^e^ihorizontal) 
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10 -11 and 10 -12 from left 



to right. The short-dashed lines represent the results obtained negleting thermal effects; the long-dashed 
lines give the baryon asymmetries if only thermal masses are included. The shaded region corresponds to 
the observed baryon asymmetryll3. 



whose flatness is lifted by the non-renormalizable interaction dTfl) and by by various super- 
gravity and finite-temperafjim corrections. Recently, leptogenesis based on this flat direction 
has been studied in detai£al!3. The complete scalar potential during an inflationary phase 
with Hubble parameter Hi is given byEj 



V(<p,T) = MT) 2 M 2 + a g a 2 s T* In l|L + ^ (^ + cc) + _1_^ 



AM 2 ' 



-ciHf\tp\- 



Hi 

8Mi 



[CLHf + C.C. 



(12) 



Here a g — 0(1) is a parameter which depends on the particle content of the theory, TO 3 / 2 a m 
reflects supersymmetry breaking at T = 0, whereas cjHf and Hjclh are due to the breaking 
of supersymmetry during inflation. Note, that the coefficient cj has to be positive for 
Affleck-Dine leptogenesis to work whereas the simplest canonical Kahler potential would 
yield the opposite sign for c/! 

The assumed phase difference between the complex numbers a m and an is responsible 
for the generation of a lepton asymmetry. During the inflationary phase the field ip is driven 
to one of the discrete minima 

\<p\ = ifo - {MHj) 1 ' 2 , arg(<p) = ^{-arg(a H ) + (2n + l)n), n = 0...3. (13) 

After inflation, when [i(T) 2 starts to be the dominant mass term, the homogeneous field <p 
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Figure 4. m„ei/e[eV] as function of sin 2 SO ; */ Am 2 oJ [eV] for the large angle MSW solution and \U e 3\ < 
0.100. 



begins coherent oscillations which eventually lead to an asymmetry in the number densi- 
ties for ordinary leptons and scalar leptons and, finally, via the sphaleron processes to an 
asymmetry in baryon number. It is remarkable that the final baryon asymmetry is rather 
insensitive to the reheating temperature after inflation in the range Tr ~ 10 5 . . . 10 12 GeV. 
A detailed calculation yieldaij, 

-3/2 fm z / 2 \a r . 



s 



t m Vl \ - 3 / 2 ( m 3 / 2 \a m \ \ 
V10- 8 eV/ V 1Te V / 



-11 x / ""vi 

6eff llO^eV 



(14) 



i.e., the baryon asymmetry is determined by the lightest neutrino mass and the strength 
of supersymmetry breaking. For an effective CP- violating angle 5 e ff ~ 0.1... 1 and 
W3/2 |<2m | — ITeV, the observed baryon asymmetry is obtained for an 'ultralight' neutrino, 

m Vl ~ (0.1 - 1) x 10~ 9 eV . (15) 

This small neutrino mass reflects the required flatness of the D-flat direction for successful 
Affleck-Dine leptogenesis. 

The prediction of an ultralight neutrino has direct implications for neutrinoless dou- 
ble beta decay. The effective neutrino mass which determines the decay rate is given by 



Y-U 2 



where m Ui are the neutrino masses and U a i the neutrino mixing 



matrix, respectively. For hierarchical neutrinos with an ultralight lowest state m Ve u e can be 
expressed in terms of the oscillation parameters determined from the solar and the atmo- 
spheric neutrino deficits, 



1 - \U e3 \ 2 ) sin 2 9 sol JAm 2 sol + \U e3 \ 2 e la J Am 2 a 



(16) 



It is remarkable that for the large angle MSW solution m„ e „ e reaches the value of m U2 ~ 
yjAm 2 ol ~ 3 x 10~ 3 — 10~ 2 [cV] (cf. fig. ||), which may be accessible in future neutrinoless 
double beta decay experiments. 



3 Thermal Leptogenesis 

Here one starts from a thermal distribution of heavy Majorana fermions which have CP 
violating decay modes into standard model leptons. The natural candidates are the three 
right-handed neutrinos Pm, i = 1 . . . 3, whose interactions are described by the lagrangian, 

C = l L (j>*X*yR-\v c R Mv R + h.c. (17) 

M is a Majorana mass matrix. The vacuum expectation value of the Higgs field (</>) = v 
generates the Dirac mass term mp = h v v which is assumed to be small compared to the 
Majorana mass M—This yields light and heavy neutrino mass eigenstates according to the 
seesaw mechanismc2l, 



with masses 



Vjv L + v c L V: , N~v R + S Ri (18) 



-Vjrn^—m D V v , m N ~ M , (19) 



where V v is the neutrino mixing matrix. 

We shall restrict our discussion to the case of hierarchical Majorana neutrino masses, 
Mi -C M2,M%. The baryon asymmetry is then determined by the CP violating decays of 
the lightest Majorana neutrino N\ = v R \ + v Rl = N, 

r(jv-ty) = i(i + e)r, r(N-+U) = ±(i-e)r. (20) 

Here T = (A^A)nAf/(87r) is the total decay width, and the parameter e <C 1 measures 
the amount of CP violation. The generation of the baryon asymmetry takes place at a 
temperature T ~ M = M\ <C M 2 , M 3 . It is therefore convenient to describe the system by 
an effective lagrangian where the two heavier neutrinos have been integrated out, 

C = IiaPKiN + N T XnCl Li (f> - ^MN T CN 

+ -jVijlLi<t>ClLj(t>+ ^ 

with 

3 



VijlIitClLrf + -vI'Im^CI^* , (21) 



^=2^ X ^Tf X kj- ( 22 ) 

fc=2 k 

can 



The CP asymmetry e arises from one-loop vertex and self-energy correctionscl(E3t3. It 
be expressed in a compact form asc-H 

3 MAV)i, Jtf _ (23) 



16tt (AtA)i 



Given the Dirac and Majorana neutrino mass matrices the CP asymmetry e and, con- 
sequently, the baryon asymmetry are determined. Consider as an example a pattern of 
fermion masses based on the group SU(5)out x U(1)f, where U(l)p is a spontaneously 



broken generation symmetry. The Yukawa couplings arise from non-renormalLzable inter- 
actions with a gauge singlet field <J? which acquires a vacuum expectation valuetj, 

Qi+Q 3 



A, 



9ij 



A 



(24) 



Here <7y are couplings 0(1) and Qi are the U(l)p charges of the various fermions, with 
Q$ = — 1. The interaction scale A is usually chosen to be very large, A > Agut- An 
example of possible charges Qi is given in table 1. The assignment of the same charge to 
the lepttp, doublets of the second and third generation leads to a neutrino mass matrix of 
the forn£§E2, 



(25) 



This structure immediately yields a large v^ — v T mixing angle. The phenomenology of 
neutrino oscillations depends on the unspecified coefficients which are 0(1). The flavour 
mixing parameter e is chosen to be ($) /A = e ~ 1/17, which corresponds to the mass ratio 

One easily verifies that the mass ratios of heavy and light Majorana neutrinos are given 

by 




Mi : M 2 : M 3 



1 



mi : m,2 '■ mz 



e 2 :l 



1 



(26) 



The masses of the two eigenstates v^ and v-j, depend on the unspecified factors 0(1). They 
are therefore consistent with the mass differences AmJ „ ~ 10~ 5 — 10 -3 eV 2 inferred from 
the large angle MSW solution of the solar neutrina,problem and Ara 2 2J , 3 ~ 10 -3 — 10 -2 eV 2 
associated with the atmospheric neutrino deficitfil. In the following we we shall use for 
numerical estimates the geometric mean of the neutrino masses of the second and third 
family, m v = [m V2 m U3 ) 1 '' 2 ~ 10 -2 eV. The choice of the charges in table 1 corresponds 
to large Yukawa couplings of the third generation. For the mass of the heaviest Majorana 
neutrino one then finds 



M, - — 



v 



10 15 GcV 



(27) 



which implies that B — L is broken at the unification scale Aqut- 

The GE asymmetry e (g3() leads to a lepton asymmetry in the course of the cosmological 
evolutionEl 



til 



Yr 



Here the factor k < 1 represents the effect of 
one has to solve the full Boltzmann equationsE 9 ''' 



.9* 



(28) 



put processes. In order to determine k 
Important processes are the AL = 2 



Tabic 1. Chiral charges of charged and neutral leptons with SU(5)gut x U(l)p symmetry. 
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Figure 5. Time evolution of the heavy neutrino number density and the lepton asymmetry for the 
SU(5)cut X U(1)f model. The solid line shows the solution of the Boltzmann equation for the right- 
handed neutrinos, while the corresponding equilibrium distribution is represented by the dashed line. The 
absolute value of the lepton asymmetry Yl is given by the dotted ^H(ie and the hatched area shows the 
lepton asymmetry corresponding to the observed baryon asymmetry 
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lepton Higgs scatterings mediated by heavy neutrinos since cancellations between on-shell 
contributions to these scatterings and contributions from neutrino decays and inverse decays 
ensure that no asymmetry is generated in thermal equilibrium. Further, due to the large 
top-quark Yukawa coupling. jane has to take into account neutrino top-quark scatterings 
mediated by Higgs bosonst3E3. These processes are of crucial importance for leptogenesis, 
since they can create a thermal population of heavy neutrinos at high temperatures T > M. 
In order to obtain the baryon asymmetry predicted by the model of neutrino masses 
described above we first have to evaluate the CP asymmmetry. Since for the Yukawa 
couplings only the powers in e are known, we will also obtain the CP asymmetry and 
the corresponding baryon asymmetry to leading order in e, i.e., up to unknown factors 
0(1). Note, that for models with a U(1)f generation symmetry the baryon asymmetry is 
'quantized', i.e., changing the U(1)f charges will change theJparyon asymmetry by powers 
of &3. One easily obtains from eqs. (|23|), (|24| ) and the tableEj, 

E ~A«.. (29) 

Using e 2 ~ 1/300 and g* ~ 100 one then obtains the baryon asymmetry, 



Y n 



K 10" 



(30) 



For k ~ 0.1 . . . 0.01 this is indeed the correct order of magnitude! The baryogenesis tern- 



perature is given by the mass of the lightest of the heavy Majorana neutrinos, 

T B ~ Mi ~ e 4 M 3 ~ 10 10 GeV . (31) 

For this model, where the CP asymmetry is determined by the mass hierarchy of light and 
heavy Majorana neutrinos, baryogenesis has been studied in detail in re£3. It is remarkable 
that the observed baryon asymmety is obtained without any fine tuning of parameters, if 
B — L is broken at the unification scale Agut- The generated baryon asymmetry does not 
depend on the flavour mixing of the light neutrinos, in particular the v^ — v T mixing angle. 

The solution of the full Boltzmann equations is shown in fig. pj. The initial condition 
at a temperature T ~ lOAfi is chosen to be a state without heavy neutrinos. The Yukawa 
interactions are sufficient to bring the heavy neutrinos into thermal equilibrium. At temper- 
atures T ~ Mi the familiar out-of-equilibrium decays set in, which leads to a non-vanishing 
baryon asymmetry. The final asymmetry agrees with the estimate ( |30| ) for k ~ 0.1. The 
dip in fig. H is due to a change of sign in the lepton asymmetry at T ~ M\. 

The final baryon asymmetry is usually obtained from the final lepton asymmetry by 
means of eq. (pi). Note, that this is conceptually not correct! The sphaleron processes are 
in equilibrium whereas the heavy neutrino decays out of equilibrium. Hence, any generated 
asymmetry is immediately distributed among other degrees of freedom in the plasma. This 
effect reduces the generated baryon asymmetry by a factor 0(1 )Ej. 

4 Towards the Theory of Leptogenesis 

The generation of a baryon asymmetry is an out-of-equilibrium process which is generally 
treated by means of Boltzmann equations. A shortcoming of this approach is that the 
Boltzmann equations are classical equations for the time evolution of phase space distribu- 
tion functions. On the contrary, the involved collision terms are S'-matrix elements which 
involve quantum interferences of different amplitudes in a crucial manner. Clearly, a full 
quantum mechanical treatment is highly desirable. This is also required in order to justify 
the use of the Boltzmann equations and to determine the size of corrections. 

All information abjaul the time evolution of a system is contained in the time dependence 
of its Green functionsL_rt3 which can be determined by means of Dyson-Schwinger equations. 
Originally these techniques were developed for non-relativistic many-body problems—More 
recently, they have—also been applied to transport phenomena in nuclear mattero, the 
electroweak plasmaEZrL3 and_th£ QCD plasmaEj. Alternatively, one may study the time 
evolution of density matriceafljLiJ. In the following we shall describe how the Green function 
technique can be used ta_qbtain a systematic perturbative expansion around a solution of 
the Boltzmann equationso. 

The time evolution of an arbitrary multi-particle lepton-Higgs system can be studied by 
means of the Green functions of lepton and Higgs fields. For the heavy Majorana neutrino 
one has 

iG a0 {x 1 ,x 2 )=Tr(^pfN a (x 1 )N (x2)) , (32) 

where T denotes the time ordering, p is the density matrix of the system, the trace extends 
over all states, and the. time coordinates t\ and t 2 he on an appropriately chosen contour C 
in the complex planec3. G(xi, x 2 ) can be written as a sum of two parts, 

G(xi,x 2 ) = @{h - t 2 )G > (x 1 ,x 2 ) + Q(t 2 - t 1 )G < {x ll x 2 ) . (33) 
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G ^H ~2iLw,22 "I • ( 34 ) 



For a system in thermal equilibrium at a temperature T = 1//3 the density matrix is 
p = exp (—pH), where H is the Hamilton operator. In this case the Green function only 
depends on the difference of coordinates and it is convenient to work with the Fourier 
transform G(p). The contour C can be chosen as a sum of two branches, C = C\ U C2, 
which lie above and below the real axis. The time coordinates are real and associated with 
one of the two branches. Correspondingly, the Green function becomes a 2 x 2 matrix, 

'G ii (p)G i2 ( P y 

G 21 {p)G 22 { P ) t 

The off-diagonal terms are given by 

G 12 (p) = G<(p), G 21 (p) = G>(p), (35) 

the diagonal terms of the matrix (p3) are the familiar causal and anti-causal Green functions. 
The free Green functions are explicitly given by 

iG>{p) = (Q(p ) - Q(p )f N (E) - 0(- P o)f N (E)) PN (p) , (36) 

iG<(p) = (0(-po) - e( Po )f N (E) - Q(- Po )fx(E)) PN (p) , (37) 
with the spectral density 

p N {p) = 2ntf + A^C-'Sip 2 - M 2 ) , (38) 
and the Fermi-Dirac distribution functions 

f N (E) = f R (E) = ^^ , E=^M 2 +p 2 . (39) 
Since N(x) is a Majorana field one has Jn — ftf, and the charge conjugation matrix C 



occurs in the spectral density (38). The Green functions S(xi,X2) and A(xi,X2) for the 
lepton doublets and the Higgs doublet also depend on chemical potentials. The Schwinger- 
Dyson equations for the Green functions G > (xi,X2), G < (xi,X2), etc. are usually referred 
to as Kadanoff-Baym equations. 

For processes where the overall time evolution is slow compared to relative motions 
the Kadanoff-Baym equations can be solved in a derivative expansion. One considers the 
Wigner transform for G(xi, X2), 

G{x,p)= /'d 4 ye i ^G(x + |,x-|) , (40) 

and for S(x,p) and A(x,p), respectively. To leading order in the derivative expansion the 
Kadanoff-Baym equations become local in the space-time coordinate x. Keeping to zeroth 
order only the on-shell part of retarded and advanced Green functions and self-energies one 
obtains the equations 

cty + j-M)G>{x,p) = cty + j-M)G<{x,p) 

- i (E> (x, p)G< (x, p) - S< (x, p)G> (x, p)) , (41) 

(~$ + !t)S>(x,p) = (~$ + lft)S<(x,p) 

= ±{ll > (x,k)S < (x,k)-n<(x,p)S > (x,p)) . (42) 
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Figure 6. One- and two-loop self energies for the lepton doublet. 



The corresponding self-energies are shown in hg. (pi). Solutions of these equations yield 
the first terms for the non-equilibrium Green functions G > (x,p)...S < (x, k) in an expansion 
involving off-shell effects and space-time variations which include 'memory effects'. 

Given the lepton self-energies we can now look for solutions of the equations ( p[ ) and 
(|42J). A straightforward calculation shows that the right-hand side of these equations van- 
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ishes for equilibrium Green functions and self-energies. Since leptogenesis is a process close 
to thermal equilibrium the solutions should be linear in the deviations, 

SG(t, P ) = G>(t,p) - G>(p) = G<(t,p) - G< (p) , (43) 

6S(t,p) = S>(t,p) - S> q ( P ) = S<(t,p) - S<(p) . (44) 

One may also expect that 8G(t,p) and 8S{t,p) can be obtained from the equilibrium Green 
functions by a small change of the distribution functions, 

iSG(x,p) = ~Sf N (x,p)p N (p) , iSS(x,k) = -e(k )5fi(x,k)pi(k) . (45) 



This ansatz indeed reduces the matrix equations (41) and ([42|) to the system of ordinary 
differential equations 



E-Sf N (t,p) = -E-f N (p) - 2(A+A) n J d$i2(p)6f N (t,p)p-pi , (46) 

d , 



gik^-Jh{t,k) = Ai m(A t 7/A * )llM / d$ l2 (k)6f N (t,pi)k- Pl 

57T 



-2(A+A) U / d<i> 1 - 2 (k)(Sf l (t,k)f 4> (p 1 ) + f l (k)5U(t,p 1 ))k-p 2 
c (AU)f/ 



M 2 



d*i53(*)(2((y/,(*,fc)^(pi) + fi(k)SU(t, Pl ) 

+8fl(t,P2)U(P3) + fl(P2)Sf ct> (t,p 3 ))k-p2 

+(W,fc)/;(Pi) + /*(fc)«,Pi) 

+28U(t,p2)U(p 3 ))k- Pl ) . (47) 

Here <?; = 6 is the number of 'internal' degrees of freedom for three generations of lep- 
ton doublets, d$i2(p) etc. denote phase space integrations, and /jv(p)> /z(p) an d /^(p) 
are the equilibrium distributions of heavy neutrinos, lepton doublets and Higgs doublet, 
respectively. 



Eqs. (|46| ) and (|47|) are the Boltzmann equations for the distributions functions of heavy 
neutrinos and lepton doublets with the correct matrix elements for decays, 2 — > 2 processes 
and CP violation. The effect of the Hubble expansion is included by means of the substi- 
tution d/dt — » d/dt — Hpd/dp. Integration over momenta then yields the more familiar 
form of the Boltzmann equations for the number densities. Note, that also the equilibrium 
distributions are time dependent since the temperature varies with time. The first term in 
eq. ( f47| ) drives the generation of a lepton asymmetry; the remaining terms tend to wash out 
an existing asymmetry. 

We conclude that a solution of the Boltzmann equations generates a solution of the 
Kadanoff-Baym equations to leading order in the expansion described above. Various cor- 
rections such as relativistiv effects, 'memory effects', off-shell effects and deviations from 
kinetic equilibrium can now be systematically studied. By means of such an analysis one 
can obtain quantitative constraints on the parameters M, (X*X)n and e, which relate the 
cosmological baryon asymmetry and neutrino properties. 
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5 Outlook 

Detailed studies of the thermodynamics of the clcctroweak interactions at high temperatures 
have shown that in the standard model and most of its extensions the electroweak transition 
is too weak to affect the cosmological baryon asymmetry. Hence, one has to search for 
baryogenesis mechanisms above the Fermi scale. 

Due to sphaleron processes baryon number and lepton number are related in the high- 
temperature symmetric phase of the standard model. As a consequence, the cosmological 
baryon asymmetry is related to neutrino properties. Generically, baryogenesis requires 
lepton number violation, which occurs in extensions of the standard model with right- 
handed neutrinos and Majorana neutrino masses. In detail the relations between B, L and 
B — L depend on all other processes taking place in the plasma, and therefore also on the 
temperature. 

Although lepton number violation is needed in order to obtain a baryon asymmetry, it 
must not be too strong since otherwise any baryon and lepton asymmetry would be washed 
out. Hence, leptogenesis leads to upper and lower bounds on the masses of the light and 
heavy Majorana neutrinos, respectively. 

Different realizations of leptogenesis imply different constraints on neutrino masses. 
Affleck-Dine leptogenesis based on the lepton-Higgs D-flat direction predicts an ultralight 
neutrino with mass m vi ~ 10~ 9 cV. Also thermal leptogenesis requires neutrino masses 
significantly below 1 eV. However, in this case the neutrino masses may be as large as 
~ 10~ 2 cV, which corresponds to the mass splitting indicated by the atmospheric neutrino 
anomaly. It is very remarkable that the observed baryon asymmetry ub/s ~ 10~ 10 is 
naturally explained by the decay of heavy Majorana neutrinos, with B — L broken at the 
unification scale Agut ~ 10 16 GeV, and in accord with present experimental indications 
for neutrino masses. 

Further work is needed to develop a full quantum mechanical description of leptogenesis 
which goes beyond the Boltzmann equations. Also important is the connection with other 
lepton and quark flavour changing processes in the context of unified theories. Finally, 
the realization of the rather large baryogenesis temperature Tg ~ I0 10 GeV in models of 
inflation should have implications for dark matter and the cosmic microwave background. 
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